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^ ■ ABSTRACT 
> ■ 

o : 

^ I Based on optical, IR and X-ray studies of Cas A, we propose a geometry 

^ I for the remnant based on a "jet-induced" scenario with significant systematic 

I departures from axial symmetry. In this model, the main jet axis is oriented 

r-| I in the direction of strong blue-shifted motion at an angle of 110 - 120° East 

^1 of North and about 40 - 50° to the East of the line of sight. Normal to this 

2 I axis would be an expanding torus as predicted by jet-induced models. In the 

to ! proposed geometry, iron-peak elements in the main jet-like flow could appear 

<^< '_ "beyond" the portions of the remnant rich in silicon by projection effects, not 

I I the effect of mixing. In the context of the proposed geometry, the displacement 

'i> I of the compact object from the kinematic center of the remnant at a position 

^ I angle of ~169° can be accommodated if the motion of the compact object is near 

^ I to, but slightly off from, the direction of the main "jet" axis by of order 30°. In 

^ J I this model, the classical NE "jet," the SW "counter-jet" and other protrusions, 

particularly the "hole" in the North, are non-asymmetric flows approximately 
in the equatorial plane, e.g., out through the perimeter of the expanding torus, 
rather than being associated with the main jet. We explore the spoke-like flow 
' in the equatorial plane in terms of Rayleigh- Taylor, Richtmyer-Meshkov and 

' Kelvin-Helmholz instabilities and illustrate these instabilities with a jet-induced 



O 



simulation. 



Subject headings: ISM: Individual: Name: Cassiopeia A, ISM: Supernova Rem- 
nants, Stars: Supernovae: General, Stars: Supernovae: Individual: Alphanu- 
meric: SN 1987A, Hydrodynamics, Instabilities 



1. Introduction 



The star we now recognize as the supernova remnant Cassiopeia A probably exploded in 
1680 (IThorstensen et al.ll200ll : iFesen et al.ll2007l ). It has been an astrophysical mystery ever 
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since its re-discovery by radio astronomers (iRyle fc Smithlll948l ) and has be en the target of 



incre a singly sophisticated studies over a range o f waveleng ths in the radio (lAnderson et al. 
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2001 



Keohane et al 



Fesen et al. 



1996 



2001 



20061). and X-ray (Mar 



Hwang et all 12000 



2006a 



Liszt fc Lucas 



199! 



Morse et al. 
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Laming et al.l l2006l ) . The projected 



19831: iHughes et all I2OO0I: IwiUingale et al.l I2OO2I . I2OO3I : 



Fesen 



Ennis et al. 



image of Gas A shows a prominent jet-like structure in the NE and a "counter-jet" diamet- 
rically opposed in the SW. Other protrusions and interruptions in the images are notable 
in images at the various wavelengths. The first image taken with the Chandra Observatory 
dramatically revealed the long-sought compact object. This dim object has yet to be given a 
firm physical identification as a neutron st ar or black hole. Th e compact object is displaced 
from the kinematic center of the remnant (IFesen et al.ll2006al ) to the South, nearly perpen- 
dicular to the locus of the "jet/counter-jet" structure. X-ray studies of the remnant have le d 
to the discovery that iron-rich regions lie beyond silicon-rich regions (IHughes et al.ll2000l ). 
suggesting some sort of inversion of the expected "onion-skin" structure of the progenitor 
star. Fine-scale turbulence might mix the originally stratified composition and lead to Fe at 
larger radii than Si, but it is not clear that such a process can lead to a large scale segregation 
of the compositions, as the observations suggest. 



Laming et al.l (120061 ) consider the possibility that the NE/SW "jets" represent a spher- 



ical explosion into an inhomog eneous GSM w i th "ca vities" so that expansion is faster into 
those cavities, as suggested by iBlondin et al. Jl996h . They find that such a model cannot 
give a sufficiently high density of plasma at high enough temperature. The ejecta expand 
so rapidly that the density is too low either for appreciable electron-ion equilibration to 
raise the electron temperature or to sufficiently ionize the plasma as obse rved for the NE 



j et. L aming et al. consider a jet-induced model based on the calculations of iKhokhlov et al. 



(119991 ) and find that such a model can provide the required conditions of ionization age and 
electron temperature for the collisional ionization equilibrium knots at the jet tip and the 
non-ionization equilibrium further back in the jet stream. Laming et al. argue that there 
should be a substantial amount of cold plasma at the "jet" head that has been cooled by 
radiative and adiabatic losses. This material could be composed of Fe, but is not neces- 
sarily so. Laming et al. note that the blast wave is not seen in the direction of the NE 
"jet." The NE jet is thus not o verdense compar e d to the stellar envelope, as it was for the 
jet-induced supernova models of [Khokhlov et al. see also lHoflich. Khokhlov fc Wang! 

2OO1I : IKhokhlov fc Hoflich|[200ll l Laming et al. find that the kinetic energy in the NE "jet" 
is a relatively modest 10^° ergs, perhaps insufficient to represent the major power source 
that caused the supernova explosion. This leaves open the basic physical mechanism of the 
explosion. 
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Hwang &: Lamingl (120031 ) argue for an eject a mas s of about 2 MfD a n d a t otal mass of 
the star at the time of the explosion of about 3 Mq. IChevalier fc Oishil (120031 ) have made 
the case for a clumpy wind that may account for the quasi-stationary flocculi and that the 
explosion mi g ht ha ve corresponded to a Type Iln or Type lib supernova, in current typology. 
Young et al.l (120061 ) argue that the progenitor was a star of 15 - 20 Mq on the main sequence 
that lost its hydrogen envelope to a binary companion and exploded with a mass of 4 - 6 

Among the elements that must be explained to fully account for the shape of Cas A are: 



A basic filametary morphology (Fig. la) 



Quasi-stationary flocculi (lFesenll200ll ) 



Fast-m oving knots of oxygen and their concentration at certain angles (IFesen et al. 



2006bf ) 



The appearance of "rings" at various parts of the perimeter (Fig. lb) 



The asymmetric distribution of mass (jWillingale et al.ll2003l ) 



The distribution of silicon and other intermediate mass elements (IHwang et al.ll2004l ) 



The concentration of iron to the SE and NW, especially its appearance at larger pro- 
i ected radii toward the SE than silico n and its apparent absence in the NE "jet" 
(iHughes et allboool : Ewang et al.lbooi l 



A prominent blue-shi fted region to the SE (iMarkert et al.l Il983l ; iHwang et al.l 12001 
Wilhngale et~aDl2002h 



The p rominent "jet" and "counter-jet" structure in the NE/SW (lFeserul200ll : lHwang et al. 



200 



The location and direction of motion of the compact object (IFesen et al.ll2006al ) 



The lack of a prominent pulsar wind nebula around the compact object (IFesen. Pavlov &: Sanwal 



2006|) 



The exp ansion velocities and e i iergetics as a func t ion of distarice, angle, and com- 



posit ion (jWillingale et al.l 120021 : iHwang et al.l l2004l : iFesen et al.l l2006al : iLaming et al. 



20061 ) 



The presence or absence of a companion star (IFesen. Pavlov fc Sanwa 
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As a step toward accounting for these properties we consider the physics and morphology 
imphed by a "jet-induced" supernova model. The motivation for such a model has arisen in 
the context of recent developments in the study of t he spectropolarirnetry of core-collapse 
supernovae flWang et al.lboosi : iMaund et all[2007al jbll3: IWang fc Wheelejboosi ) that show that 
core core collapse supernovae are routinely aspherical and frequently display a prominent 
axisymmetry. One obvious way to induc e s uch an asymrn e try is to explode the supernova 
with bi-polar jets (IKhokhlov et al.lll999l ). I Wheeler et al.l (120071 ) point out that while bi- 
polar "jet-like" flow is common in core collapse supernovae, evidence for non-axisymmetric 
structure is also prominent. The presence of non-axisymmetric morphology is manifested as 
"loops" in the plane of the Stokes parameters Q and U. This non-axisymmetric structure 
was prominent in SN 1987A (jCropper et al.lll988l ) and is now recognized to be composition- 
dependent, with different species ejected in different directions. Deeper study of Cas A may 
help to inform our interpretation of the spectropolarimetric data and vice versa. A basic 
picture of a bi-polar explosion with significant non-axisymmetric flow may apply in general 
to core-collapse explosions. Here we explore whether such a picture can help to illuminate 
the morphology of Cas A. 



2. Proposed Geometry 



As an organizing principle, we will adopt here the ansatz that the explosion of Cas A 
was jet-induced without going into detail as to the physical origin of the jet(s). A basic 
jet-induced model has certain generic features: 1) principal high-velocity flow along an axis, 
presumably the rotational axis of the progenitor star and of the new-born compact remnant; 
2) bow shocks generated by the jets, and a convergence of those bow shocks and associated 
flow onto the equator with 3) the subsequent expulsion o f the bulk of the stellar core and 



mantle in a toroidal confi guration (IKhokhlov et al.l Il999l : iHofiich. Khokhlov &: Wang 12001 



Khokhlov fc Hoflichll200ll ). For simple models in which the "up" and "down" jets are iden- 
tical, the resulting fundamental geometry comprises a jet and torus structure. Here we will 
explore what orientation of this model, and what departures from its basic depiction, may 
account for the observed structure, kinematics, and composition distribution of Cas A. A 
key component of the jet-induced m odels is the distribution of the original "onion-skin" 
layer s of the massive star pro genitor (IHofiich. Khokhlov fc Wangl l200ll : IKhokhlov fc Hofiich 
200ll : iMaeda fc Nomotdl2003l ). One expects that the jet will be predominantly composed of 
iron-peak matter arising from deep within the progenitor. The bulk of the star, composed 
of the outer layers of Si, Ca, O, He, and perhaps some H, that surrounded the iron core will 
be compressed and expelled in the expanding equatorial torus. 
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If a jet /torus geometry pertains t o Cas A, a principal question to ask is the direction of 
the main jet flow. iLaming et al.l (120061 ) have established that the the NE "jet" does represent 
a jet-like flow, but its rather feeble energy suggests that this may not be the principal axis 
of the explosion. We note that in the jet-induced model, the energy originally injected in 
a relatively narrow solid angle in the jets is redistributed throughout the stellar envelope. 
The amount of energy remaining in the original jet direction is model dependendent. Thus 
how me asurements of energy and velocities in a given direction constrain models is not 
obvious. iMarkert et al.l (119831 ) presented Einstein data that showed a substan tial blue shift 



i n the SE direction as projected on the sky, not along the classic NE "jet." iHwang et al. 



( 120011 ) and lWillingale et al.l ( 20021) presented t he sam e feature in Chandra and XMM Newton 



X-ray Doppler shift maps. Willingale et al. (|2003 ) identify a torus with an axis oriented 
approximately in this direction. Dewey et al. (j2006 ) present an other version of a Ca s A 



X-ray Doppler map coupled with transverse proper motions from iDeLaney et al.l (120041 ) to 
yield a 3D map of a selection of X-ray emitting knots. This map (comprising only 17 points, 
but promising the richness to come) shows again that, as viewed from the top, normal to 
the line of sight, the major blue-shifted component is oriented at about 45° to the East 
(clockwise) from the line of sight to the observer, with a red shifted component at about 
225°. From the observer's point of view, the knots with highest blue shift form a rough "ring" 
with a rather large opening angle of ~ 30°. This again suggests a major flow pattern at a 
projected position angle (counterclockwise from North on the plane of the sky) of about 125° 



(Hughes et al. 


2000; 


Hwang et al. 


2004) 



we assume that the direction indicated by these observations represents the projected main 
axis of the jet in the jet/torus geometry. We suggest that the main jet axis is then oriented 
to the SE at a position angle of about 125°, roughly 40 - 50° to the East and 20 - 30° to the 
South with respect to the line of sight. 

This geometry suggests an alter native explanation for the observation of iron-peak mat- 
ter beyond the silicon-rich material (iHughes et al.ll2000l ). Rather than a literal mixing of the 
matter through Rayleigh- Taylor or other processes, the apparent distribution of matter of 
different composition could be affected by projection effects. In particular, in the proposed 
geometry, the main jet would be iron rich. The equatorial flow would tend to be Si rich (also 
O and Ca rich). In this conflguration, the iron in the jet in the SE direction approaching the 
observer would tend to be faster, but it would also appear to be "outside" the Si in the SE 
portion of the slower moving torus because of projection effects. This might give a natural 
explanation for the apparent "overturn" of the Fe compared to the Si, without requiring 
"overturn," per se. 



Because of entrainment effects there could be some high- velocity Si in the axial jet flow 
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and some Fe in the equatorial flow. Note also that the initial energy in the jet will be 
dissipated throughout the mantle and envelope and that rapid flow in the direction of the 
jet could lead to adiabatic expansion and cooling so that a signiflcant portion of the jet 
flow would be difficu lt to detect in current X-ray observations, as found for the NE "jet" by 



Laming et al.l (120061 ) . 



3. The Compact Object 



Fesen et al.l (j2006al ) have shown that the position of the compact object in Cas A is 
located to the South at a position angle of 169 ± 8.4° from the kinematic center of the 
explosion determined by tracing the proper motion of expanding optical emission knots (see 
Fig. la). This has posed a special puzzle for a picture in which the classic NE jet and 
SW counter-jet represented the principal axis of the explosion. One would then expect, 
on general grounds under the jet-induced ansatz, the motion of the compact object to be 
roughly along that axis. Instead the motion inferred by the displacement of the compact 
object from the kinematic center is nearly normal to the NE/SW axis. In the geometry 
proposed here, the PA of the compact object requires some, but only modest, misalignment 
from the main axis of the jet in the SE direction. 

Consider the geometry of Cas A in spherical polar coordinates, with the observed posi- 
tion angles projected onto the plane of the sky. In spherical polar coordinates, let 9 be the 
latitudinal angle, measured from East (counterclockwise) of North and be the longitudinal 
angle measured in the plane of the observer, such that cos(0) = 1 is toward the observer and 
cos(0) = —1 is in the opposite direction with increasing to the East (clockwise as viewed 
from the North). With i,j,k unit vectors corresponding to the observer direction (0 = 0°), 
East (0 = 90°), and North, respectively, the position vectors of the jet (subscript J) and 
compact object (subscript C) are then: 

r J = sin 6 J cos + sin 6j sin 0jj + cos 0jk, (1) 

vc = sin 6c cos 0ci + sin 6c sin (pc] + cos 6c^- (2) 
The angle, x? between the jet and compact object is given by: 

rj ■ rc = |rj||rc| cosx, (3) 

such that 



cos X = sin 6 J cos j sin 6c cos (pc + sin 6j sin 0j sin 6c sin (pc + cos 6 j cos 6c. (4) 
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The observed position angle on the sky, Op^c, measured East of North (counterclockwise), is 
related to the latitudinal angle, 6c, of the 3D model of Cas A by: 



tan 6*. 



P,C 



tan 6c sin 



(5) 



Equation [5] can be solved to give sin^c^ and cos 0c in terms of sin^^ and cos6'c. When 
substituted into Eqn. HI this gives cosx as a function of cos 6'c: 

1 /2 

cos X = sin 6 J cos (pj [l — (l + tan^ 6p^c) cos^ 6c] + (sin 6j sin 0j tan 6p^c + cos 6 j) cos 6c. 

(6) 

Figure [2] gives x as a function of the unknown angle, 6c, for a range of plausible values of 6j 
and 0j for the nominal value of the position angle of the compact object, 6p^c- Note that 6c 
is bounded on one side by motion essentially on the line of sight corresponding to 6c = 90° 
for which = and cos x = sin 6j cos 0j. On the other extreme the motion at the observed 



90° and 6, 



c 



6 



p,C, 



position angle could be in the plane of the sky corresponding to (pc 
for which cosx = sin 6* j sin 0j sin + cos cos 6'p,c'- While there are ranges of parameter 
space with large values of x, there are clearly ranges that give a modest value, x ^ 30°. Fig. 
[2] shows that for a given choice of the jet angles, x has a minimum at a specific value of 
cos 6c- Figure E] gives a plot of Xmin versus 6c for 6j = 120° and a range of values of 0j and 
6p,c- 

Figures [2] and [3] show that while there are valid portions of parameter space for which 
the separation angle would be large, there is an ample and reasonable range of parameter 
space for which the direction of motion would be near to that of the proposed main jet. 
In particular, if the main axis of the jet is about 40 - 50° East and 20 - 30° South of 
the line of sight {6j = 110 — 120°), then the motion of the compact object would only 
have to be of order 30° off the main axis of the jet to point at a projected position angle 
^p,c ~ 169°. In this proposed geometry, the PA of the compact object would be nearly, if 
not exactly, aligned with the axis of the main SE jet and nearly orthogonal to the NE/SW 
"jet" axis. While this departure of the motion of the compact object from the principal axis 
requires further explanation, it is com pletely in keeping w ith the tilt of the rings in SN 1987A 
from the apparent axis of symmetry (jWang et al.l |2002| ) and with th e dispersion of prope r 
motions of young pulsars from the axes of symmetry of their remnants iNg fc Romanil (120061 ) , 
and perhaps with growing evidence from spectropolarimetry for "tilted-jet" models of core 
collapse wherein the jets de fining bi-polar axes are not exactly aligned with the dominant 
geometry of the progenitor ( Maund et al 2007a b 3). 
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4. Jets and Holes 

What, then, are we to make of the classic NE "jet" and SW "counter-jet" structure? Be- 
side these features, Cas A shows other evidence for "rings" and "holes" in optical imaging as 
spectacularly displayed in new HST ACS images as displayed in the Hubble Heritage collec- 



tion (Fig. lb; Fesen, private communication 2007; http:/ /heritage. stsci.edu/2006/30/index.html) 
In particular, there is a prominent ring in the North (on the far side from the observer) and 
a smaller ring within that, suggesting the rims of "holes." There are other smaller, but dis- 
tinct, ring- like features in the SE and SW. Clearly, Cas A does not correspond to a simple 
axially-symmetric, jet-induced flow. 

The outward flow and shocks involved in a supernova explosion are, in many circum- 
stances, subject to Rayleigh- Taylor (RT) and Richtmyer-Meshkov (RM) instabilities. In ad- 
dition, any hydrodynamic flow that is not strictly spherically symmetric will involve shear and 
the possibility of Kelvin- Helmholz (KH) instabilities. All such instabilities could contribute 



to the breakdown of strict axial symmetry. In the jet-induced models of iKhokhlov et al. 



( 119991 ) the jet was, in fact, subject to the KH instability, but the growth time was long com- 
pared to the propagation time out of the core, so the structure, although computed in full 
3D, maintained the axisymmetric structure of the initial conditions despite the presence of 
the instability. Anothe r way of breaking symmetry is to have axial jets of different strength. 



Khokhlov et al.l (119991 ) computed two axial jets of identical nature, but discussed the pos- 
sibility that the "up" and "down" jets could have different energy and momenta, with the 
difference in the momenta being delivered to the compact object in the form of a "kick." 
This imbalance in let properties could bo th break the mirror symmetry that was imposed in 
the calculations of iKhokhlov et al.l (119991 ) and cause a deviation of the directions of the jets 
and the recoil of the compact object, with possible implications for Cas A and other core 
collapse supernovae. 

We thus also hypothesize that while the main jet and counter-jet that actually triggered 
the explosion of Cas A are to the SE/NW, the classical NE/SW "jet/counter-jet" structure 
and the northern "hole" are secondary effects resulting from instabilities in the equatorial 
plane of relatively modest energy ~ 10^° ergs. In particular, we postulate that these features 
represent faster "spoke" or "flnger" "jet-like" flow in the plane of the expanding torus due 
to various instabilities encountered by the expanding toroidal flow. 

In the jet-induced model, the equatorial, toroidal component will inevitably be subject to 
KH instabilities as it propagates out into the mantle and envelope. The growth time is shorter 
for smaller wavelengths, but the shortest wavelengths will be suppressed by any restoring 
force as might be rendered by embedded magnetic flelds. It is thus difficult to determine 
the characteristic scale of these streaming instabilities in the absence of an appropriate 3D 
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numerical calculation. iHungerford et al.l (120031 ) performed relevant 3D calculations with both 
jet-like and disk-like flow, but it is not clear that their SPH calculation had the resolution 
to see the KH instability and others to be discussed below. The KH instability will also 
accompany RT and RM instabilities. There may also be a coupling between the shear flow 
along the top and bottom faces of the expanding torus and the RT and RM instabilities 
that will be triggered along the outer rim of the torus. It may be that the KH instability is 
mostly responsible for entrainment, not the growth of spokes. 

The equatorial, toroidal flow will also be subject to RT instabilities as the denser torus 
is decelerated by the less dense mantle and envelope of the star. The core is roughly an 
n = 3 polytrope. The torus expands with density decreasing roughly as if there is little 
increase in the vertical height. If the density proflle in the core is steeper than r~^, the torus 
will always tend to be denser than the mantle into which it expands. While the details will 
depend on the nature (density, energy, opening angle) of the axial jets, the calculation of 
Khokhlov et al.l (119991 ) shows that the toroidal structure remains denser than the surrounding 
He core by about a factor of 3 as the torus comes to the edge of the helium core. The torus 
will be denser than any surrounding hydrogen envelope of essentially constant density. The 
torus will thus constantly be subject to RT instabilities as it propagates outward, depending 
on its heat content and capacity to do PdV work, just as for the spherical expans i on th at 
has traditionally been studied in this context (I Chevalierl 1 1 9 76l : iKifonidis et al.ll2003l . |2006| ). 



While in the linear limit, small scales will grow faster by RT instability. In the non-linear 
limit the velocity of the growth of RT instabilities is given approximately by: 

(aL)V^ 



vrt 



(7) 



where a is the effective acceleration and L is the characteristic scale length (jYoungslll986l ). 
At this stage, larger scales will grow faster. Laser-induced production of 3D RT instabilities 
has shown that RT flngers tend to form and grow at nearly the effective "free-fall" rate, the 
velocit y of the interface before deceleration, in this case Vt, the velocity of the expanding 
torus jDrake et al.ll2004h . 



The characteristic length scale for the RT instability is often taken to be the pressure 
scale height in the star (j Che valieii 1 1 9 76l ) . Alternatively, the characteristic length scale for the 
RT structures may be the effective vertical height of the expanding torus. If so, the torus 
might tend to break up into ^ r/h fi ngers, where r is the ra dius of the torus and h its (full) 
thickness. The jet-induced models of iKhokhlov et al.l (119991 ) give an opening angle of ~ 30° 
as the torus impinges on the outer helium core. This would give r/h ~ 5 — 6, implying 5 or 6 
"fingers" per hemisphere or 10 - 12 around the full perimeter of the torus. This is somewhat 
larger than, but perhaps in the ball park of, the number of "fingers" seen in Cas A. The 
most prominent fingers might be fewer than the total. 
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Equation [7] depends on the amplitude of the deceleration, which is itself a function 
of the density of the torus and the ambient medium into which it propagates. Invoking 
momentum conservation (neglecting the pressure within the torus as if it were already cold 
and expanding homologously) and expressing the effective acceleration as the ram pressure 
exerted by the mantle, ~ PmanVt divided by the mass per unit area of the torus that is 
decelerated in a given time, one can write: 



a ~ 



PtVtT 



Pman'^t 

Ptr 



(8) 



where pman is the density of the mantle, pt is the density of the torus, Vt is the velocity of 
the torus, and r is the timescale of the growth of the RT structure. From Eqns. [7] and [HI we 
can write 

Pman'^t 



T 



pta 



1/2 



and hence 



givmg 



pt 



2 2 

L' 



Vrt 



pt 

where L ~ r^/S and is the radius of the torus. 



-Vt, 



(9) 



(10) 



(11) 



For the particular simulation of the core of a 15 M© star of iKhokhlov et al.l (119991 ). the 
torus hits the edge of the helium core of radius 2 x 10^° cm in 35 seconds and hence has a 
mean velocity of ~ 6000 km s~^. The torus has a density about 3 times that of the mantle 
so the acceleration that drives the RT instability at that point will be a ~ 10'' cm s~^ and 
the growth rate of the fingers will be vrt ~ I'^t- The latter result suggests that the RT 
fingers might not have reached the saturation limit at this epoch . Note that if the fingers 
grow at the rate of motion of the interface, Vt ( jPrake et al.l 120041 ). then they could develop 
~ 3 times faster. 

The RM instability occurs when shocks propagate down density gradients and especially 
when shocks encounter density discontinuities, such as found at composition boundaries in 
massive stellar cores. For the cores of stars that were the likely progenitor of Cas A, these are 
the Fe/ Si boundary, the Si/0 bou ndary, the C-O/He boundary, and the He/H boundary, 
if any (IKifonidis et al.l l2003l . |2006| ). For perturbations of initially small amplitude, 6r, of 
characteristic length, L, the growth rate for RM instabilities is: 



Vrm = ^TTVt 



6r 



A, 



(12) 
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where A is the Atwood number, A 



pt pman ( Youngsi 1 19861 ) ■ At the boundary of the hehum 

pt~rpman ' 

core, where A ~ 1/2 then vrm ~ '^vt ~ dv^T for 6r ~ L. 



To illustrate the capacity of a jet-induced model to produce this variety of instabihties 
(if not literally the structure of Cas A), we show in Figure H] the result of a two-dimensional 
simulation of equal and opposite jets propagating from the center of a helium star. The 
progenitor model wa s the 3.5 Mq core remaining from the evolution of a star of 10 Mq 
(IWoosley et al.lll995l ). The jets were introduced as an inflow boundary condition at the 
inner radial coordinate of 3.82x10^ cm with an opening angle of about 25°. The jet density, 
pressure and injection velocity rise linearly to maximum values of 6.5x10^ g cm~'^, 10^^ ergs 
cm~^ and 3.22x10^ cm s~^, respectively, over the course of 0.05 seconds. After about 0.5 
seconds, the jet injection velocity is linearly reduced to zero in about 1.5 seconds. These 
parameters were chosen so as to deliver a total energy of 10^^ ergs to the star (one-half 
that energy in each jet). The mass interior to the lower boundary, 1.593M0, was included 
as a point mass source of gravitation. The self-gravity of the remaining portion of the 
star was neglec t ed. T hese conditions were intended to approximately reproduce those of 



Khokhlov et al.l t9m . 



The resulting dynamics were computed using the FLASH code (IFryxell et al.l l2000l ) . 
We used a grid in spherical coordinates with eight levels of reflnement in both radius and 
angle. The maximum effective number of angular zones is 1024 and the maximum number of 
radial zones is 8192. This gives a minimum resolution in radius of ~ 1.5 x lO'' cm and in angle 
of 0.1 75° (corresponding t o 3 x 10^ cm at R = 1 x 10^ cm). This compares to the simulation 
of iKhokhlov et al.l (119991 ) that had a minimum resolution of ~ 3.7 x 10^ cm near the inner 
boundary that degraded in the Cartesian grid of that simulation to ~ 2.3 x 10^ cm at the 
outer boundary (~ 1.5 x 10^^). Our maximum resolution remains constant i n radius and ou r 



angula r resolution is higher throughout the simulation compared to that of iKhokhlov et al. 
(I1999I ). The instabilities that set in at larger radius are thus better resolved. 

Figure m shows the results after 47 seconds with the jets in the vertical directions. The 
leading shock induced by the jets has propagated out of the star (initial radius 1.07 x 10^^ 
cm) and off the grid (1.2 x 10^^ cm). The effect of the jets in the inner core have "healed" 
due to transverse pressure gradients. The remnants of the passage of the jets can still be 
seen in the outer triangular regions at a radius of 4 — 7 x 10^ cm. The effects of instabilities 
are plainly seen. The RM instabilities and associated KH instabilities induced as the jet 
shock propagated across the Si/0 boundary are seen in the "cap" at about 9 x 10^° cm. The 
horizontal structure is in the equatorial plane and results from the convergence of the jet- 
induced blast waves on the equator. KH rippling is seen at about 5 x 10^° cm and an RM/KH 
"mushroom" is seen on the leading edge of the toroidal structure at about 8 x 10^° cm. All this 
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structure would be different for 3D (the small extensions right on the computational axis at 
9 X 10^° cm are surely artificial) and with a true, rather than numerically generated, Reynold's 
number in the flow, but the basic presence of the instabilities is well illustrated. We note that 
while a plausible, but arbitrary, perturbat ion must be added to induce RM and RT growth 
in a spherical model (IKifonidis et al.l 120061 ). a jet-induced model, by assumption, represents 
a large asymmetric "perturbation," leading to "natural" and significant instabilities. Details 
of this simulation will be presented elsewhere. 

The RT and RM (and KH for that matter) instabilities tend to lead to mixing in the 
non- linear limit. If we want to explain "fingers" in Cas A, then it may be necessary for the 
fingers to break out of the star before that turbulent, non-linear limit is reached, although 
finer scale structure might be related to the fast-moving knots. The analysis given above 
(Eqns 8-11) suggests that might be the case for the RT instability; that the RT fingers have 
not reached saturation before the explosion leaves the helium core. The RM instabilities 
may grow faster locally at the composition discontinuities, but it is not clear in the absence 
of an appropriate simulation how this will interact with the RT structure. Limitation of the 
non-linear turbulence might be aided by the fact that the progenitor star had lost most or 
all of its hydrogen envelope and hence was something similar to a bare helium core. Figure 
m which is not yet in homologous expansion, hints that some of this structure may survive. 

We note that in the proposed picture the toroidal flow and "fingers" should be composed 
of elements such as Si, Ca, O, and should be relatively devoi d of iron-peak matt er. This 
might agree with the high- velocity oxygen-rich knot structure of iFesen et al.l (j2006bl ) and the 
relative paucity of Fe in the direction of the NE "jet" observed in the X-ray. Once again, we 
note that especially fast moving matter in the "fingers" may undergo adiabatic expansion 
and cooling so that the bulk of it may be hard to detect. 

If there are equatorial "fingers" of faster flow, this flow may also drive lateral pressure 
waves piling up matter on the "walls." These walls might thus be denser and that higher 
density might inhibit the flow adjacent to the primar y flow in the finger. This might account 
for the marked gaps in the flow pattern reported by lFesen et al.l (l2006bl ). 



Another way to induce RT instabilties is to accelerate denser material by lighter mate- 
rial. This might occur if there is a later injection of energy from the new-born neutron star 
into the previously exploded and expanding material. There are a variety of mechanisms 
that might provide a somewhat delayed input of energy into increasingly less dense matter 
surrounding a new-born neutron star. 



Wheeler &: Akiyamal (120071 ) have pointed out that non-axisymmetric instabilities in 



rapidly rotating neutron stars are likely to generate a strong magnetoacoustic flux. This 
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flux would push on the inside of the expanding torus produced by a previous phase of MHD 
jets that pr opagated up the rota t ion a xis. The resulting interaction might produce RT in- 
stabilities. IWheeler fc Akiyamal (120071 ) note that this process is likely to be amplified by 
the deleptonization, contraction and spin-up of the proto-neutron star on the timescale on 
which the radius contracts, of order 0.1 - 1 s. It is not clear how any such magnetoacoustic 
flux would be propagated: some might go up the rotation axis, but a substantial portion 
might go out the equatorial plane, directly accelerating the postulated expanding torus from 
within. 

At somewhat later phases, after the neutron star has shrunk to its final radius, ~ 
10 km, it will continue to cool by neutrino emission for another 10 - 100s, the (other) 
Kelvi n- Helmholtz timescale. During this phase the neutron star may emit a neutrino-driven 
wind (iQian fc Woosleylll996l ). Toward the end of this phase, as the density declines in the 
vicinity of the neutron star , the flow may become magn etically-do minated, thus resulting 



in a relativistic win d flow (IBucciantini et al.l |2006| . 120071 ) (see also [Thompson et al.l 12004 



Metzger et al.ll2007l ). IBucciantini et al.l (120071 ) note that this sort of flow seems to naturally 
occur in the equatorial plane. This work does not account for the fact that the supernova 
ejecta will not have dispersed on the relevant timescales and thus will present a different 
outer bound ary condition - a working surfa ce - than that assumed in the numerical models as 
discussed by iKomissarov fc Barkovl (120071 ). Nevertheless, there is a suggestion for a natural 
tendency for an original MHD jet flow up the axis to evolve to an equatorial flow at later 
times driven by a pulsar wind. The characteristic timescale, 10 - 100 s is roughly the time 
it takes for the original shock, 10 s, and the torus, somewhat slower, to propagate to the 
edge of a helium core. Thus, at about the time this jet/torus structure breaks out of the 
core, there might be a strong, relativistic, equatorial pulsar wind that begins to accelerate 
the jet-induced torus from within. The low densities associated with this late-time fast wind 
will naturally lead to RT instabilities as the slower, denser torus is accelerated from within 
by the pulsar wind. We argue that such a wind should be expelled into structure roughly 
similar to Figure HI 



In effectively free-streaming calculations, IBucciantini et al.l (120071 ) flnd that their nu- 
merical models give essentially the same asymptotic solutions as analytic solutions with 
the energy scaling like sin^ 6. This would mean that half the energy is directed within a 
ha lf- angle of ~ 50 °. Thi s energy distribution is wider than the torii computed here and 



by iKhokhlov et al.l (119991 ). but an appropriate calculation must consider the interior of the 
supernova as the appr opriate "outer" boundary condition. This might lead to a more con- 
flned equatorial flow. IBucciantini et al.l (120071 ) also note a tendency to form "plasmoids" 
with characteristic angles of 5 - 10°. Such structures r night see RT instabilitie s and give 
somewhere of order 20 - 40 "flngers." The calculations of IBucciantini et al.l (120071 ) were done 
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for a magnetar-like case and the results for the density and Lorentz factor of the pulsar wind 
will depend on the specific magnetization and inner boundary condition, but the principles 
should apply to both pulsars and magnetars. We note that arguments have been made that 
the compact objec t in Cas A is a magnet ar, but also that there is no c urrent obvious evidence 
for a pulsar wind (IFesen. Pavlov fc Sanwalll2006l : iKrause et al.l 120051 ) . 



5. Discussion and Conclusions 

Despite attempts to deconvolve the three-dimensional structure, the interpretation of 
the kinematics, dynamics, and morphology of Cas A remains difficult. We adopt a jet- 
induced structure as an organizing principle to frame the analysis of observations and to 
provide a model to test. We suggest that Cas A exploded from a jet-like mechanism, but 
that, unlike many interpretations, the main jet is not the well-known NE "jet." We argue that 
a structure to the SE with a position angle of ~ 125° represents the main jet. This feature is 
less immediately discernable because of projection effects, but is suggested in Doppler maps. 
It is also iron rich. In the jet-induced model, the bulk of the ejecta, the intermediate-mass 
elements, should be ejected in an equatorial torus (for jets of equal and opposite momentum) 
that results from the convergence of blast waves on the equator. We interpret the NE "jet" 
and other structure as flows that are roughly normal to the main SE jet that result from 
non-axisymmetric structure in the equatorial plane of the jet. This structure may resolve 
one of the current mysteries of Cas A, that the proper motion of the compact object seems 
to be nearly normal to the principle axis of the dynamics. In our proposed morphology, the 
compact object has been "kicked" at an angle ~ 30° to that of the main jet. The proposed 
picture also opens the way for a new interpretation of chemical structure, especially the fact 
that iron is seen beyond silicon: this observation may be significantly affected by projection 
effects, rather than resulting from true "mixing" of the ejecta. 

There have been h ints of the geometry we propose here in previous discussions of Cas 



A. iMarkert et al.l (Il983l ) proposed an inclined ring model (their Figure 5), but their solution 
was a ring with its axis virtually along the line of sight. To account for the SE feature, 
they invoked non-uniform emissivity around the perimeter of the ring. The ring may have 
non-uniform emissivity, but their particular solution does not account for why this SE fea- 
ture is iron-rich, and does not seem to be quantitatively consistent with subsequent work. 



Willingale et al.l (120021 ) used XMM Newton composition and Doppler data to do a 3D decon- 



volution, presenting a proposed side view from the East (their Fig ure 11). They remark that 



the data is well characterized by the doughnut shape suggested by lMarkert et al.l (119831 ) . but 



their side view seems to show the axis of the torus tilted at about 45° to the line of sight. 
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rather than the 87° ad vocated by iMarkert et aLl (119831 ) . With the spatial resolution avail- 



able, the project i on of IWillingale et al.l (|2002| ) shows no particular evidence of the NE jet 
Willingale et al.l (120031 ) propose a clumpy toroidal geometry in which the axis of the toroid 
is oriented 40° East of the line of sight and 125° (counterclockwise) to the South of North, 
essentially the same as we have chosen here. Along this torus they identify concentrations 
of mass and energy on opposing sides to the "North" and "South" that are near to, but 
not quite the same as North and South on the plane of the sky (see their Figure 4 for the 
coordinate system they define). They refer to these concentrations of ma ss and energy as 



"jets" but do not relate them to the Doppler feature in the SE discussed by IWillingale et al. 



( 120021 ). and they are distinctly not the directions of the classic NE and SW jets. It is als o 
not quite clear how th is torus relates to the deconvolution given by IWillingale et al.l (120021 ). 
Willingale et al.l (120031 ) suggest that the torus they identify may be the shell of matter ejected 
by a progenitor Wolf-Rayet star. Here we raise the possibilit y that the torus i s tha. t pro- 
duced by a jet-induced explosion. The mass map presented by IWillingale et al.l (120031 ) does 
not show any special evidence for the iron-rich mass in the SE that we identify as the main 
direction of the jet axis. This may be rela ted to the fact that th e hot iron moving along 
the jet axis has relatively little mass. Like IWillingale et al.l ( 120031 ) . we have no immediate 
explanation for the asymmetric distribution of mass around the perimeter of this torus, but 
propose that this is related to the instabilities and non-axisy mmetric flows normal to the jet 
axis that we have discussed here, rather than "jets" per se. IWillingale et al.l (120031 ) point 



out that the apparent asymmetry of the explosion might suggest shear and hence turbulence 
and dumpiness intrinsic to the ejecta without needing a collision with an external medium. 
We have illustrated just this sort of effect in Figure HI 



Krause et al.l (120051 ) have provided another interesting and mysterious aspect concerning 
the nature of Cas A. They have identified an apparent bi-polar flow and associated infrared 
echoes. Their preferred orientation is a position angle of about 26° at an angle of about 
82° to the line of sight, nearly in the pla ne of the sky. This angle differs from every other 
geometric feature discussed in this paper. iKrause et al. J2005h suggest that Cas A contains 
a magnetar that had a soft gamma-ray repeater outburst circa 1952. There is no sign of a 
Crab-like synchrotron nebula, which is also an issue for our suggestion that a pulsar wind 
might play a role in shaping the ejecta. The energy injected by a soft gamma-ray repeater 
burst is rather small compared to a supernova kinetic energy, but such a burst of energy 
might affect portions of the ejecta, enhancing previous irregularities. It would be interesting 
to look for other possible indications of s uch a magnetar burst. We note that the ionization 
ages identified by IWillingale et al.l (120031 ) are on the scale of a century or less. 



Burrows et al.l (120051 ) proposed that the main axis of the explosion of Cas A is that 



associated with the "iron-rich, mass-rich, energy-rich" SE direction and that this direction 
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be a rotation axis. They speculate that the rotation axis precessed from the generation 
of the main jet to the generation of the later, secondary flow associated with the NE jet, 
or that the latter could be a magnetic dipole axis (although a magnetic axis would be 
expec ted to rotate in direc tion at the pulsar period if tilted with respect to the rotation 
axis). iBurrows et al.l (120051 ) argue that the inferred direction of motion of the central object 
is along the proposed axis of explosion with an inferred blue-shift. We have made the speciflc 
proposal of a displacment of the line of motion of the compact object by ~30° from the axis of 
the main jet. iBurrows et al.l (120051 ) suggest that the m otion of the c ompa ct object could be 
due to a slightly top-bottom asymmetric bi-polar flow. IJanka et al.l (120051 ) also suggest that 
the NE jet is induced after the main explosion, perhaps due to accretion onto the neutron 
star. Here we propose that the re is an explici t jet /t orus structure of t o Cas A, and agree 
with the general suggestions of IBurrows et al.l (120051 ) and iJanka et al.l (120051 ) that the NE 
jet is some form of secondary flow, caused by the complex interaction of the jets or by the 
late-time effect of a pulsar wind, or perhaps repeated soft gamma-ray repeater outbursts. 

To produce a robust explosion, the main explosion engine should involve ~ 10^^ ergs. 
Even if this energy begins in jets, it is spread through a large solid angle by subsequent 
dynamical interaction with the surrounding progenitor star. Thus it is not clear how much 
energy should be directed along the main axis once free expansion is attained. This can be 
quantifled by suitable r nodels. It rem a ins to be seen whether an analyis of the SE structure 
similar to that done by iLaming et al.l (120061 ) will reveal a sufficiently large kinetic energy in 
the SE to be consistent with a jet in that direction being the primary origin of the explosion. 
The evidence that only ~ 10^° ergs is involved in the NE jet is consistent with that being 
only a secondary effect. It is also important to determine the extent to which projection 
of an iron-rich SE jet on a silicon-rich equatorial flow can quantitatively account for the 
observed apparent composition inversion. 

An important aspect of this model is that the non-axisymmetric structure arises from 
instabilities in or near the equatorial plane. We have sketched possible instabilities, including 
the effects of a puls ar wind. Three-dimensio nal hydrodynamics calculations or even pehaps 
MHD calculations ( IStone fc Gardined 120071 ) are required to determine whether or not the 
structure of Cas A can arise in a natural and physically self-consistent way from a jet-induced 
model. Such a model would have to yield a pla usible explanation o f transverse velocities up 
to 14,000 km s"^ in the direction of the NE jet (IFesen et al.ll2006bl ) in contrast to transverse 
velocities of ~ 8,000 km s~^ along the SE structure (the radial velocities in both of these 
directions are more modest; IWillingale et al.l (120021 )). We note that this high velocity in 
the NE direction may be an indication that this is a secondary flow, perhaps driven by a 
pulsar wind, invoking once again the notion that the original jet energy must be dissipated 
throughout the ejecta to cause the explosion. A separate injection of energy may be necessary 
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to blow holes in t he original e j ecta, a ccelerate some regions to high velocity and produce 



fast moving knots. iFesen et al.l (l2006bl ) note that the NE jet is like a "spray" of ejecta knots 



rather than a narrow jet. That might be an important clue to the underlying process. 

There is no model today that accounts for the nature of Cas A. The positive aspects of 
the model we present here are that it roughly accounts for the iron-rich SE structure, the 
extension of iron beyond silicon in that direction, and the motion of the compact object in 
nearly that direction. On the negative side, there is much structure in Cas A that remains 
unexplained. Our model gives no natural explanation for why the NE jet and SE counter-jet 
seem so nicely and oppositely aligned, a problem for any model in which these are secondary 
flows. It may also be the case that the structure seen in Cas A might be generated in the 
absence of a jet-induced origin. This has certainly not been demonstrated. It may also 
be true that the NE jet defines a symmetry axis and the compact object has been kicked 
at some oblique angle due to equatorial instabilities, but such a theory would also require 
elaboration. 

An important aspect going forward is to compare models with observations. For our 
model we predict: a) there is enough mass and kinetic energy in the SE and NW directions 
to account for the explosion (some of this matter may be fast, cool and difficult to directly 
observe; much of the originally directed energy has spread throughout the ejecta); b) there 
is high- velocity iron extending to large distances in the SE (in front) and NW (in back); 
c) the 3D structure and kinematics should show an expanding toroidal structure, roughly 
normal to the SE direction we have defined here, consistent with the energy deposited in 
the progenitor star by the jet-induced blast waves; d) the intermediate mass elements are 
concentrated in this torus; e) the apparent "overturn" that places iron beyond silicon in the 
plane of the observer is largely due to differentiated, composition-dependent, bulk flow, not 
mixing due to turbulence; f) the compact object is moving toward us; g) the "holes" in the 
remnant caused by secondary flow, including the NE and SW jets, are the locus of faster 
moving material; h) the walls deflning the holes are in transverse motion; the holes should 
be getting larger. We urge those analyzing the data on Cas A to put these predictions to 
the test in order to more deeply unders tand this remark a ble ex plosion. The structure to the 



SE (and NW) and the torus deflned by IWillingale et al.l (120031 ) . in particular, warrant much 



more detailed study. We need to better understand the kinematic and dynamical relation 
between the SE feature, the NE/SW "jets," and this torus and to know whether the torus 
is more consistent with a pre-existing shell or a structure created by the supernova. 
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Fig. 1. — (a - left) X-ray image of Cas A and illustration of the coordinate system for this 
work. The solid lines represent the observed position angle of the compact object, Op^ci 
and that of the proposed main jet, Op^j, with dashed hnes indicating a plausible range of 
angles for the jet. Heavy dashed lines illustrate the approximate location of the NE and 
SW "jet" and "counter-jet" structure that we interpret here as a secondary flow resulting 
from instabilities. Adapted from Hwang et al. (2004). (b - right) Optical image of Cas A 
illustrating some of the "holes" that characterize the morphology of Cas A. Adapted from 
http: / /heritage.stsci.edu/2006/30/index.html 
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cos{6^) 

Fig. 2. — The angle, x, between the axis of the proposed main jet and that of the motion 
of the compact object is given as a function of the direction of recoil of the compact object, 
9c, for plausible choices of the orientation of the jet axis, 9j and (pj. Angles 9 are measured 
East (counterclockwise) from North and angles are measured East from the observer line 
of sight (clockwise as viewed from the North). Horizontal lines mark given values of the 
separation angle, x — 10, 20, 30 and 40°. The range of cos 9c is restricted to fall between 
and -1, since the compact object is known to be recoiling toward the SE and that is also the 
assumed direction of the main jet in this work. The angle 9c is bounded on one side by 9c 
= 90° for which (pc = and cos x = sin 9j cos 0j and on the other extreme by (pc — 90° for 
which 9c — 9p^c and cos x = sin 9j sin 0j sin 9p^c + cos 9j cos 9pfi. 
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Fig. 3. — The minimum value of the separation angle, Xmm, between the axis of the proposed 
main jet and that of the motion of the compact object as shown in Fig. [2] is given for the 
orientation of the jet axis 9j = 120° and a range of plausible values of (pj, and for the 
nominal observed value of the position angle of the compact object, 9p^c = 169°, and 1, 2, 
and 3 a limits on that angle where a = 8.4°. Angles 6 are measured from the North and 
angles (p are measured from the observer line of sight. Horizontal lines mark given values of 
the separation angle, x = 10, 20, 30 and 40°. 
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Fig. 4. — (a - left) The result of launching sjTumetric vertical jets with 10^^ ergs into a helium 
star of 3.5 Mq. This 2D simulation was performed with the FLASH code with a central point 
source of gravity, but neglecting the self-gravity of the remaining matter. The computation 
was done on a 180° grid and a mirror image used to create this figure. The light strip down 
the middle is an artifact of this reflection process. This image corresponds to 47 seconds 
after the initiation of the explosion after the leading shock has left the star. The asymmetry 
of energy injection naturally leads to a variety of hydrodynamic instabilities associated with 
the direct propagation of the jets and the secondary, equatorial toroidal flow resulting from 
the collision of the jet-induced blast waves on the equator, (b - right) Expanded view of 
the equatorial torus showing the KH "cap" and associated wrinkling on the top and bottom 
surface of the flow. See the electronic edition of the Journal for a color version of this figure. 



